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Abstract 

The steady boundary layer flow of nanofluid over an exponential stretching surface is investigated analytically. The 
transport equations include the effects of Brownian motion parameter and thermophoresis parameter. The highly 
nonlinear coupled partial differential equations are simplified with the help of suitable similarity transformations. 
The reduced equations are then solved analytically with the help of homotopy analysis method (HAM). The 
convergence of HAM solutions are obtained by plotting /i-curve. The expressions for velocity, temperature and 
nanoparticle volume fraction are computed for some values of the parameters namely, suction injection parameter 
a, Lewis number Le, the Brownian motion parameter Nb and thermophoresis parameter Nt. 
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1 Introduction 

During the last many years, the study of boundary layer 
flow and heat transfer over a stretching surface has 
achieved a lot of success because of its large number of 
applications in industry and technology. Few of these 
applications are materials manufactured by polymer 
extrusion, drawing of copper wires, continuous stretch- 
ing of plastic films, artificial fibers, hot rolling, wire 
drawing, glass fiber, metal extrusion and metal spinning 
etc. After the pioneering work by Sakiadis [1], a large 
amount of literature is available on boundary layer flow 
of Newtonian and non-Newtonian fluids over linear and 
nonlinear stretching surfaces [2-10]. However, only a 
limited attention has been paid to the study of exponen- 
tial stretching surface. Mention may be made to the 
works of Magyari and Keller [11], Sanjayanand and 
Khan [12], Khan and Sanjayanand [13], Bidin and Nazar 
[14] and Nadeem et al. [15,16]. 

More recently, the study of convective heat transfer in 
nanofluids has achieved great success in various indus- 
trial processes. A large number of experimental and the- 
oretical studies have been carried out by numerous 
researchers on thermal conductivity of nanofluids 
[17-22]. The theory of nanofluids has presented several 
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fundamental properties with the large enhancement in 
thermal conductivity as compared to the base fluid [23]. 

In this study, we have discussed the boundary layer 
flow of nanofluid over an exponentially stretching sur- 
face with suction and injection. To the best of our 
knowledge, the nanofluid over an exponentially stretch- 
ing surface has not been discussed so far. However, the 
present paper is only a theoretical idea, which is not 
checked experimentally. The governing highly nonlinear 
partial differential equation of motion, energy and nano- 
particle volume fraction has been simplified by using 
suitable similarity transformations and then solved ana- 
lytically with the help of HAM [24-39]. The convergence 
of HAM solution has been discussed by plotting h- 
curve. The effects of pertinent parameters of nanofluid 
have been discussed through graphs. 

2 Formulation of the problem 

Consider the steady two-dimensional flow of an incom- 
pressible nanofluid over an exponentially stretching sur- 
face. We are considering Cartesian coordinate system in 
such a way that #-axis is taken along the stretching sur- 
face in the direction of the motion and y-axis is normal 
to it. The plate is stretched in the ^-direction with a 
velocity U w = U 0 exp (x/l). defined at y = 0. The flow 
and heat transfer characteristics under the boundary 
layer approximations are governed by the following 
equations 
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(4) 



where 



where (u, v) are the velocity components in (x, y) 
directions, pf is the fluid density of base fluid, v is the 
kinematic viscosity, T is the temperature, C is the nano- 
particle volume fraction, (pc) p is the effective heat capa- 
city of nanoparticles, (pc)f is the heat capacity of the 
fluid, a = kl{pc)f is the thermal diffusivity of the fluid, 
D B is the Brownian diffusion coefficient and D T is the 
thermophoretic diffusion coefficient. 

The corresponding boundary conditions for the flow 
problem are 



u — U w (x) = U 0 exp (x/l) , v=-p(x), T = T WI C = C w at y = 0, 
u = 0, T = C = Coo as y -> oo, 



(5) 



in which U 0 is the reference velocity, f}(x) is the suc- 
tion and injection velocity when P(x) > 0 and fi(x) < 0, 
respectively, T w and are the temperatures of the 
sheet and the ambient fluid, C w , are the nanoparti- 
cles volume fraction of the plate and the fluid, 
respectively. 

We are interested in similarity solution of the above 
boundary value problem; therefore, we introduce the fol- 
lowing similarity transformations 



u = U 0 exp (y)f (rj) , v = -^J^-exp (^) {/ (rj) + rjf (rj)} , 



T-To, 



(6) 



Making use of transformations (6), Eq. (1) is identi- 
cally satisfied and Equations (2) -(4) take the form 



jr\r\r\ + ffr)r) 2/ - 0, 



+ Pr (fOr, - f v 0 + NWrfgr, + NtO 2 ) = 0 



Nt 

gr,r 1 +Le(fg ri -f ri g) + —O riri =0, 



f = —v W/ U = 1/ 0 = 1, g = 1 at rj = 0, 
f v — ► 0, 0 — ► 0, g — ► 0 as rj — ► oo, 



(7) 
(8) 

(9) 
(10) 



(PC)p 



Nt = D B —^ (C w - C^) , Nb - 
(pc) f 



D T (PC) P (T w - T^) . v D v 
, he = — , Pr = — . 

Too (pc)f v D B a 



The physical quantities of interest in this problem are 
the local skin-friction coefficient Cf, Nusselt number 
Nu x and the local Sherwood number Sh x , which are 
defined as 



f* = — -u> 

pU 2 Q ei 



(T w - T x ) dy I 



(C w - Coo) dy | 0 



V2^c /t = f" (0) , nu x /JirT x = -Me' (0) , sh x /JmT x = -^'(o), 
where Re^ = U^lv is the local Renolds number. 



(id 



3 Solution by homotopy analysis method 

For HAM solutions, the initial guesses and the linear 
operators L t (i = 1 - 3) are 

fo fo) = 1 - v w - e"", 0o (17) = e~\ g 0 (r,) = eW) 

A (f) =/'"-/', C2(.e) = e"-e, c 3 (g)=g"m 

The operators satisfy the following properties 
Ci [ae-* +c 2 e v +c 3 ] = 0, (14) 



C 2 [cie- 11 +c 5 e ri ] =0, 
C 3 [c 6 e-v + c 7 e^]=0, 



(15) 



(16) 



in which C x to C 7 are constants. From Equations (7) 
to (9), we can define the following zeroth-order defor- 
mation problems 

(l - p) A [/ (v, p) -fo (i;)] = ph.H.N, [f (r,,p)] ,(17) 
(1 - p) £ 2 p (r/,p) - 0 O (i;)] = p^H 2 N 2 [§ (77,p)] (18) 
(1 - p) C 3 [g (77, p) - go (77)] = ph 3 H 3 N3 [g (17, p)] , (19) 
f(0,p) = -v w , f(0,p) = l, f(oo,p) = 0, (20) 



0(O,p) = l, <9(co,p) = 0, 
g'(0,p) = l, g(oo,p) = 0. 



(21) 



(22) 



In Equations (17)-(22), hi, h 2 , and h 3 denote the non- 
zero auxiliary parameters, Hi, H 2 and H 3 are the non- 
zero auxiliary function (Hi = H 2 = H 3 = 1) and 
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(24) 
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2 ' d n dr,* ■ Nb" m J" Nbdrf (25) 



Obviously 

f(ri,0)=f 0 (ri), f (r/, 1) = fin), 

§(r ) ,0)=9 0 (r l ), 0 fa, 1) = 0 fa) , 

£ fa, 0)= So fa), £fa, l)=£fa). 



(26) 
(27) 
(28) 



When p varies from 0 to 1, then f (rj,p), 0 (rj,p), 
g(f],p) vary from initial guesses f 0 (77), 0 O (77) and g 0 (77) 
to the final solutions/ (77), 6 (77) and g (77), respectively. 
Considering that the auxiliary parameters hi, h 2 and h 3 
are so properly chosen that the Taylor series of f (rj, p), 
0 (rj, p) and g (rj, p) expanded with respect to an embed- 
ding parameter converge at p = 1, hence Equations (17)- 
(19) become 



/fa,p)=/ofa) + X>fa)p m , 
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00 
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ml 
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The mth-order problems are defined as follow 

A [fm fa) - Xmfm-1 fa)] = ^l^m fa) , 



(29) 



(30) 



(31) 



(32) 



(33) 



(34) 



(35) 



C 2 [0m fa) ~ Xm0m-1 fa)] = h 2 R 2 m fa) , 
£3 [gm fa) ~ Xmgm-l fa)] = Mm fa) - 

/ m (o)=/;(o)=/;(oo) = o, 

0m (0) = 0 m (00) = 0, 
g' m (0) = gm (00) = 0, 

where 



(36) 
(37) 
(38) 
(39) 
(40) 

(41) 



JO, m < 1, 
Xm " j 1, m > 1. 

m— 1 m— 1 

*m W =/»-! (^) + £/m-l-*/Z' ~ ^ E/m-l-X (42) 



fe=0 



*m fo) = Cl + ft" E {fm-1-*^ " + NK_l-kg» + Nt^-l-k^} , (43) 



iVt 

#m fo) = £-1 + ^ E {/m-l-fegfc - fm-l-M + ■ (44) 

k=0 

Employing MATHEMATICA, Equations (35)-(40) 
have the following solutions 
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(45) 



(46) 



-(47) 



in which a° 0 , n , n , n are the constants and the 
numerical data of above solutions are shown through 
graphs in the following section. 

4 Results and discussion 

The numerical data of the solutions (45)-(47), which is 
obtained with the help of Mathematica, have been dis- 
cussed through graphs. The convergence of the series 
solutions strongly depends on the values of non-zero 
auxiliary parameters h t (i = 1, 2, 3, hi = h 2 = h 3 ), which 
can adjust and control the convergence of the solutions. 
Therefore, for the convergence of the solution, the h- 
curves is plotted for velocity field in Figure 1. We have 
found the convergence region of velocity for different 
values of suction injection parameter v w . It is seen that 
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Figure 1 ft-Curve for velocity. 



with the increase in suction parameter v w , the conver- 
gence region become smaller and smaller. Almost simi- 
lar kind of convergence regions appear for temperature 
and nanoparticle volume fraction, which are not shown 
here. The non-dimensional velocity/ against r\ for var- 
ious values of suction injection parameter is shown in 
Figure 2. It is observed that velocity field increases with 
the increase in v w . Moreover, the suction causes the 
reduction of the boundary layer. The temperature field 
6 for different values of Prandtle number Pr, Brownian 
parameter Nb, Lewis number Le and thermophoresis 
parameter Nt is shown in Figures 3, 4, 5 and 6. In Fig- 
ure 3, the temperature is plotted for different values of 
Pr. It is observed that with the increase in Pr, there is a 
very slight change in temperature; however, for very 
large Pr, the solutions seem to be unstable, which are 
not shown here. The variation of Nb on 6 is shown in 
Figure 4. It is depicted that with the increase in Nb, the 
temperature profile increases. There is a minimal change 
in 0 with the increase in Le (see Figure 5). The results 
remain unchanged for very large values of Le. The 
effects of Nt on 0 are seen in Figure 6. It is seen that 



Figure 3 Variation of temperature for different values of Pr 
when Le = 2, h = -0.1, Nt = Nb = 0.5, v w = 1. 




Figure 4 Variation of temperature for different values of Nb 
when Le = 2,h = -0.1, Nt = 0.5, v w = 1, Pr = 2. 



temperature profile increases with the increase in Nt; 
however, the thermal boundary layer thickness reduces. 
The nanoparticle volume fraction g for different values 
of Pr, Nb, Nt and Le is plotted in Figures 7, 8, 9 and 10. 
It is observed from Figure 7 that with the increase in 
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Figure 7 Variation of nanoparticle fraction g for different 
values of Nb when Le = 2, h = -0.1 , Nt = 0.5, v w = 1, Pr = 2. 



Pr= 1 
Pr = 2 



Figure 8 Variation of nanoparticle fraction g for different 
values of Pr when Le = 2, h = -0.1, Nt = 0.5, v w = *\, Nb = 0.5. 




= 1 

= 1.5 



Figure 10 Variation of nanoparticle fraction g for different 
values of Nt when Le = 2, h = -OA, Nt = 0.5, v w = 1, Pr = 2. 



Nb, g decreases and boundary layer for g also decreases. 
The effects of Pr on g are minimal. (See Figure 8). The 
effects of Le on g are shown in Figure 9. It is observed 
that g decreases as well as layer thickness reduces with 
the increase in Le. However, with the increase in Nt, g 
increases and layer thickness reduces (See Figure 10). 
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